Origin of the Three-body Parameter Universality in Efimov Physics 
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In recent years extensive theoretical and experimental studies of universal few-body physics have 
led to advances in our understanding of universal Efimov physics. Whereas theory had been the 
driving force behind our understanding of Efimov physics for decades, recent experiments have 
contributed an unexpected discovery. Specifically, measurements have found that the so-called three- 
body parameter determining several properties of the system is universal, even though fundamental 
assumptions in the theory of the Efimov effect suggest that it should be a variable property that 
depends on the precise details of the short-range two- and three-body interactions. The present 
Letter resolves this apparent contradiction by elucidating previously unanticipated implications of 
the two-body interactions. Our study shows that the three-body parameter universality emerges 
because a universal effective barrier in the three-body potentials prevents the three particles from 
simultaneously getting close together. Our results also show limitations on this universality, as it is 
more likely to occur for neutral atoms and less likely to extend to light nuclei. 
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In the early 70's, Vitaly Efimov predicted a strikingly 
counterintuitive quantum phenomenon today known 
as Efimov effect: in three-body systems for which the 
two-body s-wave scattering length a is much larger than 
the characteristic range ro of the two-body interaction, 
an infinite number of three-body bound states can be 
formed even when the short-range two-body interactions 
are too weak to bind a two-body state (a < 0). The 
Efimov effect, once considered a mysterious and esoteric 
effect, is today a reality that many experiments in ultra- 
cold quantum gases have successfully observed and con- 
tinued to explore 0413 ■ One of the most fundamental 
assumptions underlying our theoretical understanding of 
this peculiar effect is that the weakly bound three-body 
energy spectrum, and other low-energy three-body scat- 
tering obscrvablcs, should depend on a three-body pa- 
rameter that encapsulates all details of the interactions 
at short distances [la ]. For this reason, the three-body 
parameter has been viewed as nonuniversal since its value 
for any specific system would depend on the precise de- 
tails of the underlying two- and three-body interactions 

EMI. 

In nuclear physics, this picture seems be consistent, 
i.e., three-body weakly bound state properties seem to 
be sensitive to the nature of the two- and three-body 
short-range interactions More recently, however, 

Berninger et al. [|[ have directly explored this issue 
for alkali atoms whose scattering lengths are magneti- 
cally tuned near different Fano-Feshbach resonances . 
Even though the short-range physics can be expected to 
vary from one resonance to another, Efimov resonances 
were found for values of the magnetic field at which 
Q=Qtk=— 9-l(2)r v dWi where r v( jw is the van der Waals 
length [2(| HI]]. Therefore, in each of these cases, the 
three-body parameter was approximately the same, thus 
challenging a fundamental assumption of the universal 
theory. Even more striking is the observation that the 
Efimov resonance positions obtained for 39 K 7 Li 0- 



7], 6 Li and 85 Rb [Hj] are also consistent with val- 

ues of a3 b /?*vdw found for 133 Cs [j|. (Note that the 
work in Ref. Q also provided early suggestive evidence 
of such universal behavior.) These observations provide 
strong evidence that the three-body parameter has uni- 
versal character for spherically-symmetric neutral atoms, 
and therefore suggest that something else beyond the uni- 
versal theory needs to be understood. 

In this Letter, we identify precisely the physics be- 
yond the universal theory needed to explain a universal 
three-body parameter, presenting theoretical evidence to 
support the recent experimental observations. Previous 
work has shown that the three-body parameter can be 
universal — that is, independent of the details of the 
interactions — in three-polar-molecule systems [22j ] and 
in three-atom systems near narrow Fano-Feshbach reso- 
nance [H, [24| , although recent work has shown that the 
latter case likely requires even more finely-tuned con- 
ditions [251 ] . Our present numerical analysis, however, 
adds another, broader class of systems with a universal 
three-body parameter: systems with two-body interac- 
tions that efficiently suppress the probability to find any 
pair of particles separated by less than ro (see Section A 
in Ref. [IB])- This class of systems, therefore, is more 
closely related to systems near broad Fano-Feshbach res- 
onances fl9| . 

Such a suppression could derive from the usual classical 
suppression of the probability for two particles to exist 
between r and r + dr in regions of high local velocity 
HkL{r), which is proportional to [hki,{r) / mdr\~ l (m be- 
ing the particle mass) . the time spent in that interval dr 
(see Section B in Ref. [26[ ) . It is possible that there could 
be an additional suppression as well, through quantum 
reflection from a potential cliff [27] ] . Systems supporting 
many bound states, such as the neutral atoms used in 
ultracold experiments with their strong van der Waals 
attraction, clearly exhibit this suppression. In general, 
a finite-range two-body potential that supports many 
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FIG. 1: (a) The full energy landscape for the three-body potentials at a = oo for our v° model potential supporting 25 
two-body bound states, (b) effective diabatic potentials W v relevant for Efimov physics for v\ with an increasingly large 
number of bound states (A^ is the value of A that produces a = oo and n s-wave bound states). W v converge to a universal 
potential displaying the repulsive barrier at R ~ 2r v dw that prevents particles access to short distances, (c)-(e) demonstrate 
the suppression of the wave function inside the potential well through the channel functions & v (R;9,ip) for R fixed near the 
minima of the Efimov potentials in (b). (c) shows the mapping of the geometrical configurations onto the hyperangles 9 and 
ip. (d) and (e) show the channel functions, where the "distance" from the origin determines |$„| 1/ ' 2 , for two distinct cases: in 
(d) when there is a substantial probability to find two particles inside the potential well (defined by the region containing the 
gray disks) and in (e) with a reduced probability — see also our discussion in Fig. [2] In (d) and (e), we used the potentials 
Wsch and Vx, respectively, both with n = 3. Note that we have used 9, (p) = <Ev(-R; tt — 9,ip) to generate the lower half of 

the surface plots. 



bound states decreases steeply with decreasing interpar- 
ticle distance r, starting when r/r v( jw 1, at which 
point the potential cliff plays a role analogous to a repul- 
sive potential for low-energy scattering. We demonstrate 
this fact by showing that the three-body parameter in 
the presence of many two-body bound states roughly co- 
incides with that for a 100% reflective two-body model 
potential, where the two-body short-range potential well 
is replaced by a hard-sphere. 

The starting point for our investigation of the univer- 
sality of the three-body parameter is the adiabatic hy- 
perspherical representation (l8l . |2&| . This representation 
offers a simple and conceptually clear description by re- 
ducing the problem to the solution of the "hyperradial" 
Schrodinger equation : 



h 2 d 2 
Y^IlR 2 



W V (R) 



F V {R) 



+ W VV ,{R)F V ,{R) = EF V {R). (1) 

Here, the hyperradius R describes the overall size of 
the system; v is the channel index; /i = m/y/3 is the 
three-body reduced mass for particle masses m; E is the 
total energy; and F v is the hyperradial wave function. 
The nonadiabatic couplings W uv > drive inelastic transi- 
tions, and the effective hyperradial potentials W v support 



bound and resonant states. To treat problems with deep 
two-body interactions — necessary to see strong inside- 
the-well suppression — requires us to solve Eq. ([1} for 
two-body model interactions that support many bound 
states, a challenge for most theoretical approaches. Us- 
ing our recently developed methodology 23], however, 
we were able to treat systems with up to 100 two-body 
ro-vibrational bound states and solve Eq. ^ beyond the 
adiabatic approximation. Here, we explore the universal- 
ity of the three-body parameter by analyzing a number 
of model potentials [see v sc h, v b x , v^dw m ^qs. (S.l)- 
(S.4) of Ref. (26[] with different numbers of bound states. 

Figure Q] (a) shows the adiabatic potentials £/„ at 
| cz, I = oo obtained using only the two-body Lennard- Jones 
potential, v^, supporting 25 dimer bound states. At first 
glance, it is difficult to identify any universal properties 
of these potentials. Efimov physics, however, occurs at 
a very small energy scale near the three-body breakup 
threshold. Indeed, a closer analysis of the energy range 
\E\ < h 2 /mrl dw [Fig. 0] (b)] reveals the universal prop- 
erties of the key potential curve controlling Efimov states 
and universal scattering properties. 

Figure [T{b) shows one of our most important pieces of 
theoretical evidence for universality of the three-body pa- 
rameter: the effective adiabatic potentials W v obtained 
using for more and more two-body bound states con- 
verge to a single universal curve. [In some cases in 
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(a) A = A; (b) A= A; 




FIG. 2: Density plot of the three-body probability density 
|<E>„ (R; 9, ifi)\ 2 sin 29 which determines the three particle con- 
figuration [see Fig. [T] (c)] in the 9-tp hyperangular plane for 
a fixed R (sin 29 is the volume element), (a)-(d) show the 
results for an R near the minimum of the Efimov potentials 
in Fig. [T] (b) for the first four poles of the v% model as indi- 
cated, (a) shows that there is a negligible probability to find 
the particles at distances smaller than r v dw (outer dashed 
circle) and, of course, inside the 1/r 12 repulsive barrier (in- 
ner solid circle). For higher poles, i.e., as the strength of 
the hard-core part of Vx potential decreases, the potential 
becomes deeper and penetration into the region r < r v dw 
is now classically allowed. Nevertheless, (b)-(e) show that 
inside-the-well suppression still efficiently prevents the par- 
ticles from being found at distances r < r v dw- In fact, we 
calculated the probability to find the atoms at r < r v< jw and 
found to be in the range 2%-4%. 

Fig. Q] (b) we have manually diabatizcd W v near sharp 
avoided crossings in order to improve the visualization 
(see Section D in Ref. [26|).] As one would expect, 
the usual Efimov behavior for the effective potentials, 
W v =-tf{s% + l/4)/2^/J? 2 with s ~ 1.00624, is recov- 
ered for R > 10r v dw- It is remarkable, however, that W v 
also converge to a universal potential for R < 10r v( jw 
and, more importantly, these effective potentials display 
a repulsive wall or barrier at R rs 2r v dw- This barrier 
prevents the close collisions that would probe the short- 
range three-body physics, including three-body forces 
known to be important in chemistry, thus making the 
three-body parameter universal as we will confirm be- 
low. This is in fact our most striking result: a sharp 
cliff or attraction in the two-body interactions produces a 
strongly repulsive universal barrier in the effective three- 
body interaction potential. 

Qualitatively, this universality derives from the re- 
duced probability to find particles inside the attractive 
two-body potential well. This effect can be seen in terms 
of the channel functions <f>„ (see Section C in Ref. 26}]), 
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FIG. 3: The Efimov potential obtained from the different 
two-body potential models used here. The reasonably good 
agreement between the results obtained using models sup- 
porting many bound states (v BC h, and Vx) and wj," w [ob- 
tained by replacing the deep potential well with a hard wall 
but having only one (zero-energy) bound state] supports our 
conclusion that the inside-the-well suppression of the wave 
function is the main physical mechanism behind the univer- 
sality of the three-body effective potentials. The differences 
between these potentials are shown to lead to differences of 
a few percent in the three-body parameter. Note that the 
two-body potentials V\, v\, and «vlw ah have an asymptotic 
van der Waals tail, whereas v BC h decays exponentially at large 
distances. 



in Figs. [T] (c)-(e) and the hyperangular probability den- 
sities in Fig. [5] In the adiabatic hyperspherical repre- 
sentation, the space forbidden to the particles fills an 
increasingly larger portion of the hyperangular volume 
as R decreases. This evolution can be visualized as the 
dashed lines in Fig. [5] (a)-(d) expanding outward. In 
the process, the channel function $„ is squeezed into an 
increasingly smaller volume, driving its kinetic energy 
higher and producing the repulsive barrier in the uni- 
versal Efimov potential. Moreover, this suppression im- 
plies that the details of the interaction should be largely 
unimportant. Consequently, different two-body model 
potentials should give similar three-body potentials. Fig- 
ure[3]demonstrates this universality by comparing W v ob- 
tained from different potential models supporting many 
bound states. Perhaps more importantly, it compares 
them with the results obtained from the two-body model 
^vdw th 8 ^ replaces the deep well by a hard wall, essen- 
tially eliminating the probability of observing any pair 
of atoms at short distances. Quantitatively, however, the 
fact that the barrier occurs only at R ss 2r v dW indicates 
that universality might not be as robust as in the cases 
discussed in Refs. |22h25| . It is thus important to quan- 
tify the value of the three-body parameter to assess the 
size of nonuniversal effects. 

In principle, the three-body parameter could be de- 
fined in terms of any observable related to the Efimov 
physics [r| . Two of its possible definitions are [l5| : the 
value of a = a^7 h < at which the first Efimov resonance 



4 




FIG. 4: Values for the three-body parameter (a) k* and (b) a^ h as functions of the number n of two-body s-wave bound 
states for each of the potential model studied here. The values for a^ h were obtained by determining the value of a at which 
the Efimov state becomes unbound, (c) Experimental values for a^ b for 133 Cs 0] (red: x, +, □, and *), 39 K [J (magenta: 
A), 7 Li (blue: •) and [3,01 (green: ■ and o), 6 Li @, @| (cyan: A and V) and 0, (brown: T and 0), and 85 Rb [lj] 
(black: ♦). The gray region specifies a band where there is a ±15% deviation from the «vdw results. The inset of (a) shows the 
suppression parameter £p n [Eq. (S.5) in Ref. [Hj]] which can be roughly understood as the degree of sensitivity to nonuniversal 
corrections. Since £p n is always finite — even in the large n limit — nonuniversal effects associated with the details of the 
short-range interactions can still play an important role. One example is the large deviation in «« found for the w sc h (n = 6) 
model, caused by a weakly bound g-wave state. For n > 10 we expect k* and a^ h to lie within the range of ±15% stablished 
for n < 10. 



appears in three-body recombination (see for instance 
Ref. [H) and = (m^l/ft 2 ) 1 / 2 , where E° h is the 
energy of the lowest Efimov state at \a\ — > oo. Our nu- 
merical results for k* and a^ h are shown in Figs.|4ja) and 
(b), respectively, demonstrating their universality in the 
limit of many bound states. In fact, the values for re* and 
a^ h in this limit differ by no more than 15% from the uJj W 
results — k* = 0.226(2)/r vc iw and a^ h = — 9.73(3)r v dw 
[solid black line in Fig. HJa) and (b)] — indicating, once 
again, that the universality of the three-body parameter 
is dependent upon the suppression of the probability den- 
sity within the two-body potential wells. We note that 
our results for single channel two-body models should 
be applied for broad Feshbach resonances. Given this 
picture, we attribute the non-monotonic behavior of re* 
and a^ h in Fig. 2] to the small, but finite, probability to 
reach short distances, thus introducing nonuniversal ef- 
fects related to the details of two- and three-body forces 
including interactions with an isolated perturbing chan- 
nel. Nevertheless, our results for Og b are consistent with 
the experimentally measured value for 133 Cs @, Hj], 39 K 
0, 7 Li [EH3, 6Li IHH IHi, and 85 Rb 0, all of which 
lie within about 15% of the v^ w result. Curiously, if one 
simply averages the experimental values, then the result 
differs from the v^ w result by less than 3%. 

Previous treatments have failed to predict the univer- 
sality of the three-body parameter for various reasons. 
In treatments using zero-range interactions, for instance, 
the three-body parameter enters as a free parameter to 
cure the Thomas collapse [32| , preventing any statement 
about its universality. Finite range models, like those 
used in some of our own treatments [Tij [corresponding 
to the results for w sc h with n = 2 and 3 in Figs. @] (a) and 



(b)], have failed for lack of substantial suppression of the 
probability density in the two-body wells. This scenario, 
however, should reflect better the situation for light nu- 
clei having few bound states and shallow attraction. In 
contrast to Ref. [18| , other models [24], [33j-[38j have found 
better agreement with experiments. Our analysis of these 
treatments, however, indicates that the two-body mod- 
els used have many of the characteristics of our w^ w , 
therefore satisfying the prerequisite fora universal three- 
body parameter. A recent attempt [39| to explain the 
universality of the three-body parameter avoided explicit 
two-body models altogether, using instead an ad hoc hy- 
perradial potential that bore little resemblance to our 
numerical potentials in Fig. [TJ This ad hoc three-body 
potential displayed strong attraction at short distances 
in contrast to our key finding, that a cliff of attraction 
for two bodies produces a universal repulsive barrier in 
the three-body system. Consequently, even thou gh a uni- 
versal three-body parameter was found in Ref. [39(, the 
fundamental understanding provided by the approach is 
uncertain (see our discussion in Section D of Ref. |26l|). 

In summary, our theoretical examination shows that 
the three-body parameter controlling many of the uni- 
versal properties of Efimov physics can be a universal 
parameter under certain circumstances that should be re- 
alized in most ultracold neutral atom experiments. Pro- 
vided the underlying two-body short-range interaction 
supports a large number of bound states, or it has some 
other property leading to the suppression of the wave 
function at short distances, three-body properties asso- 
ciated with Efimov physics can be expected to be univer- 
sal. While these arguments suggest universality also for 
the three-body parameter in heteronuclear systems that 
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exhibit Efimov physics with only resonant interspecies 
two-body interactions, verifying this prediction is a high 
priority for future theory and experiment. Equally im- 
portant is the exploration of the relationship between 
a < and a > Efimov features — currently a sub- 
ject of a number of controversies 0] — under the new 
perpective our present work offers. 
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SUPPLEMENTARY MATERIAL 

A. Two-body model interactions 

Our theoretical model for the two-body interactions 
mimics the tunability of the interatomic interactions 
via Fano-Feshbach resonances [l[ by directly altering 
the strength of the interparticle interactions and, conse- 
quently, leading to the desired changes in a. In this work, 
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we have considered various model interactions in order to 
test the universality of our three-body calculations — we 
consider a result is universal if it is independent of the 
particular choice of the underlying two-body interaction. 
The models we have used for the two-body interactions 
are: 



v sc h{r) = -Dsech (r/r ) ■ 
«S(r) = -^(l-A"/r») 
^(r)=-^exp(-A 6 /r 6 

..hs 



«E w (r) = B hs S (r hs - r) - -£0 (r - n 



(S.l) 
(S.2) 

(S.3) 

(S.4) 



The potential model in Eq. (IS.ll) is the modified Poschl- 
Teller potential, where D determines the potential depth; 
Eq. (|S.2p is the usual Lennard- Jones potential; and 
Eq. (|S.3[) is a dispersion potential with a soft wall at short 
range. In Eq. (|S.4[) . 0(x) the step-function [Q(x) = for 
x < and 1 elsewhere. In practice, however, we have 
used a smooth representation of the step-function in 
order to simplify our numerical calculations.] The po- 
tential in Eq. (|S.4[) . therefore, consists of a hard-sphere 
potential for r < rh s (-Bhs ^ ^W^L) and a long-range 
dispersion — 1/r 6 potential for r > rh s - In the present 
study, the parameters D and A in Eqs. (|S.1|) - (IS.3I) are ad- 
justed to give the desired a and number of bound states. 
For convenience, we denote the values of D and A at 
which there exist zero-energy bound states (|a| — > oo) as 
£)* and A* , where n corresponds to the number of s-wave 
bound states. For the potential model in Eq. (|S.4[) . how- 
ever, we adjusted rh s to produce the changes in a, but 
we only performed three-body calculations near the first 
pole at rh s ~ 0.8828r v dw- 

The potential models in Eqs. (IS.2I) - (|S.4[) have in com- 
mon the same large r behavior given by the van der 
Waals interaction, ~Ce/r 6 , with Cq the van der Waals 
dispersion coefficient. For these potentials, the impor- 
tant length scale is the van der Waals length defined 
as r vd w = (2/i2bC' 6 //i 2 ) 1 / 4 /2 Q, where /i 2 b is two-body 
reduced mass. Therefore, in order to compare the re- 
sults from these models to those for v sc h, we define an 
equivalent r v dw for w sc h through the relationship be- 
tween r v dw and the effective range r c g for w sc h, namely, 
r v dw ~ r e g/2. 78947 .2], valid for values \a\ ^> r . In 
fact, for u sc h we have found that r e ff(|a| = oo) depends 
on the potential depth £>*, while for van der Waals type 
of interactions, such as those in Eqs. (|S.2[) - (IS.4I) . r c g is 
fixed by r cff « 2.78947(2 / Lt 2b C 6 /ft 2 ) 1 / 4 /2. Nevertheless, 
after rescaling the three-body results obtained from v sc h 
in terms of r v dw 5 we found good agreement with those 
obtained using v b x , and f^w tnus indicating that the 
effective range is in fact the most relevant length scale in 
the problem. We chose, however, to keep our results in 
terms of r v dw to make a closer analogy to the experi- 
ments using alkali atoms. 



B. Suppression of inside-the-well probability 

Our claim here is that the origin of the universality 
of the three-body parameter is related to the suppres- 
sion of the probability to find two particles at distances 
r < r v dw • This suppression leads to the formation of the 
universal three-body potential barrier near R « 2r v dWj 
thus strongly suppressing the three-body wave function 
at shorter distances where the details of the two- and 
three-body interactions are important [see Fig. 1 (b) of 
our main text]. Next we explore the origin of this sup- 
pression at the two-body level. 

To gain some insight into the likelihood of finding par- 
ticles inside the potential well, we start by defining the 
following quantities, 

= - / I^W| 2 dr, (S.5) 

C p ut (k) = lim — L_ f \Mr)\ 2 dr, (S.6) 

where V'fcM is the two-body scattering wave function at 
energy E^b = k 2 /2/i2b, defined such that 



r^oo sin(fcr + S) 
sin S 



(S.7) 



with S = S(k) the s-wave scattering phase shift. This 
definition for V'fcM? therefore, leads to a zero-energy 
(k —> 0) wave function of the form: ipo(r) = 1 — r/a. 
[Note that in the above equations ro is the characteris- 
tic range of the two-body interaction. For the potential 
model in Eq. (|S.1[) . ro is just the quantity in the argu- 
ment of the sech function, while for the potential models 
in Eqs. (jS.2[) - (|S.4[) it is defined to be r = r vd w]- 

The parameters ^ n and £° ut can be associated with 
the "average" amplitude of the wave function inside and 
outside the potential well, respectively. We, therefore, 
define the amplitude inside the well relative to the am- 
plitude outside the well as: 



£ P out (fc) 



= 2&{k) sin 2 6. 



(S.8) 



The relative amplitude above vanishes in the limit k — > 
(sin 8 ~ —ka), as a result of our choice for the asymptotic 
solution in Eq. (|S. T[) . except at |et| = oo, when S is an odd 
multiple of tt/2. The quantity ^(k), however, remains 
finite in the k — > limit. We believe, therefore, that 
£"(£; — > 0) to be the most relevant parameter for our 
analysis on the inside-the-well suppression. Rigorously 
speaking, ^ n is not a probability, but it does measure the 
likelihood of finding two particles within r < ro, where 
the short-range interactions are experienced. Figure [STTI 
shows a typical result for |a|/ro and £™(£; — > 0) for the 
two-body potential in Eq. (|S.1|) . Figure [STTI shows that 
in the universal regime near the poles in a, the wavefunc- 
tion is suppressed (small ^ n ) and that this suppression 
becomes more efficient as the potential becomes deeper 
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and more states are bound. The black filled circles, open 
circles and open squares in Fig. IS. 11 showing the values of 
£™ at \a\ — > oo, a = 5ro, and a = — 5tq, respectively, illus- 
trate this trend. Note, however, that for values \a\ < ro, 
the parameter £™ quickly increases, indicating a higher 
likelihood to find particles inside the potential well. Sim- 
ilar results are also obtained for the potentials and v\ 
[Eqs. (fST2|) and (jST3|) . respectively]. 




4 6 8 10 12 

{mrlD/h 2 ) 1 / 2 

FIG. S.l: The red solid line represents the scattering length, 
a/ ro, while the green dashed lines represents the parameter 
Cp{k — > 0). Both quantities are given as a function of the 
depth D of the two-body interaction model v sc h [Eq. (|S.1|) ]. 
whose values for which \a\ = oo are indicated in the figure 
as D^, where n is the number of s-wave states. Black circles, 
open circles and open squares are the values of ^ n at \a\ — > oo, 
a = 5ro, and a = — 5ro, showing the suppression of the ^ 
as the number of bound states increases. The results for £p n 
also show higher efficiency of the inside-the-well suppression 
for |o|/ro 3> 1. The black dash-dotted line shows the semi- 
classical suppression factor l/fc_L(ro) demonstrating the origin 
of the suppression mechanism. 

The suppression in £™ (for a fixed value of \a\ ^> r ) 
can be understood to be a result of the usual semiclassical 
suppression of the wave function. In the WKB approach 
[3H5|, the wave function inside the potential well can be 
written as, 



?AwkbO <r ) = 



C 



(r) 



- sin 



k L (r')dr' + - 



, (S.9) 



fc|(r) 



where C is a normalization constant and 
2/i2b[-E2b — v ( r )] defines the local wave number fci(r). 
Therefore, for deep potentials, the suppression of the 
wave function inside the potential well is simply related 
to factor fci(r)" 1 / 2 which leads to the suppression of the 
amplitude of the WKB wave function [Eq. (|S.9|) ] between 
r and r+dr. Physically, this can be seen as the increase of 
the local velocity ftfcz, (r) / m (m being the particle mass) 
and, therefore, the decrease of the time spent in that 
interval dr, [hk^r) / mdr]' 1 . Therefore, in the WKB ap- 
proximation, one would expect £™ to be simply propor- 



tional to I/Zcl- We test this expectation by plotting in 
Fig. IS . 1 1 the value of l/fe/,(ro) (black dash-dotted line) 
which shows that the suppression in — > 0) is consis- 
tent with the semiclassical suppression described above. 

We have also considered the possibility that two-body 
quantum reflection [6, 7] could be responsible for the sup- 
pression of inside-the-well probability, and consequently 
the universality of the three-body parameter. Quantum 
reflection is a fundamental phenomenon involving reflec- 
tion of a low energy quantum matter wave from a poten- 
tial cliff. In the context of ultracold collisions of ground 
state neutral atoms, quantum reflection can occur when 
two atoms approach within a distance comparable to the 
range of the interaction, ro, and it is caused by the rapid 
change of the potential in this region [7j. In that case, 
just a fraction of the incoming wave probes the details 
of the atom-atom interaction at r < rg. Nevertheless, 
our analysis leads us to believe that this is of secondary 
importance. The main reason is that in a three-body 
system, even at zero (three-body) energy, two particles 
can have energies that are not low enough to ensure that 
quantum reflection plays a rule. In fact, accordingly to 
our estimates, in a three-body system the average two- 
body energy will be proportional to 1/2/ii? 2 . This im- 
plies that near the universal barrier (i? s» 2r vc iw) the 
two-body energy can be of the order of l/mr^ dW . Based 
on a model similar to the one used in Ref. [7| the reflec- 
tion probability at such energies is about 10%, which is 
too low to explain our observations. Our above analysis, 
however, shows that for such energies £p n (fc) does not dif- 
fers substantially from £™(0)i offering a more consistent 
picture than quantum reflection. Evidently, to precisely 
quantify the importance of quantum reflection requires 
a much more detailed analysis of how the energy of two 
particles in a three-body system behaves as a function of 
R, which is beyond the scope of our present study. 



C. Three-body hyperspherical adiabatic 
representation 

As mentioned in our main text, our analysis of the uni- 
versality of the three-body parameter was based on the 
hyperspherical adiabatic representation. The following, 
is a sketch of its basic aspects and fundamental equa- 
tions. A more detailed description can be found, for in- 
stance, in Refs. [8|,|9(. We start here from the hyperradial 
Schrodinger equation : 



H 2 d 2 
2p'dB 2 



W V {R) 



F^R) 

. V ,{R)F V .{R) = EF V {R). (S.10) 



The hyperradius R can be expressed in terms of inter- 
particle distances for a system of three equal masses m, 
R = 3~ 1//4 (r 2 2 + 7-33 + f'li) 1 / 2 , and describes the overall 
size of the system; v is a collective index that represents 
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all quantum numbers necessary to label each channel; 
/i = m/V3 is the three-body reduced mass; E is the to- 
tal energy; and F v is the hyperradial wave function. The 
simple picture resulting from this representation origi- 
nates from the fact that nonadiabatic couplings W uv i are 
the quantities that drive inelastic transitions, and that 
the effective hyperradial potentials W Vl like any usual 
potential, support bound and resonant states, but now 
for the three particle system. 

To obtain W v and W uv > , however, one needs to diago- 
nalize the adiabatic Hamiltonian H^d which includes the 
hyperangular kinetic energy and all of the interactions 
[H, Q ■ Our present study has assumed that the interac- 
tions are given as a pairwise sum of two-body model po- 
tentials and used those from Eqs. (|S.l|) - (|S.4p . The eigen- 
values of i? a d, determined for fixed values of R by us- 
ing the method in Ref. § , are the three-body potentials 
U V (R) [main text Fig. 1 (a)] associated to the respec- 
tive channel functions (eigenfunctions) 3>„(.R;f2) [Tol. TlTjj . 
Therefore, the channel functions <fr v contain all the in- 
formation necessary to determine, for fixed values of R, 
the geometric structure of the three particle system. For 
our present study concerning three identical bosons in 
the total angular momentum state J™ = + (it being the 
parity) , the set of hyperangles ft [9[ is reduced to simply 
{8, ip}. This helps to visualize the geometric structure of 
the three-body system [see Figs. 1 (c)-(e) and Fig. 2 of 
our main text]. 

Finally, to determine the effective potentials W v and 
nonadiabatic couplings W vv i , which are the main quan- 
tities in Eq. (IS.10|) . one can use the expressions, 



W V [R) = U V (R) - —Q VV (R), 
2/i 



W VV ,{R) = - — 



ft 2 
2/z 



2P VV ,(R)— + Q VV ,(R) 
dR 



(S.ll) 
,(S.12) 



where 



Pw{R) = Q)| ^|<M#;^)}), (S.13) 

Q VV .{R) = n)\^\^(R;n)}}. (S.14) 



In the above expressions, the double-bracket matrix ele- 
ments indicate that integrations are carried out over only 
the hyperangles J7. As is mentioned in the main text, to 
treat problems with deep two-body potentials, and con- 
sequently many bound states, we have used a modified 
version of the above formulation since P uv < and Q U v'i 
defined above, can become numerically unstable. This 
problem can be overcome using the formulation devel- 
oped in Ref. [l2j]. Apart from this more technical aspect, 
the above formulation gives the general ideas concerning 
the hyperspherical representation. We encourage the in- 
terested reader to look for more details in Refs. |8l. fol. [l2l| . 



D. Effective adiabatic potentials 

The effective potentials shown in Fig. 1 (a) of the main 
text are very complicated, making identification of the 
important physics challenging. For that reason, this sec- 
tion presents some details of our work that not only give 
support to our physical interpretation of the nature of the 
three-body parameter but also show how some important 
nonuniversal aspects appear in the hyperspherical adia- 
batic representation. In fact, Fig. IS.2I shows some of the 
most drastic nonadiabatic effects found in our calcula- 
tions. The model proposed in Ref. is also considered, 
and we show that this model offers an interesting qualita- 
tive picture; but since incorporation of corrections to that 
model diminishes its accuracy in the three-body parame- 
ter universality, its agreement with experiment might be 
fortuitous. 

We first consider the validity of the single- channel adi- 
abatic hyperspherical approximation and point out the 
manner in which some important nonuniversal features 
manifest themselves. Figure IS.2I shows the results for 
U V {R) and W V (R) obtained from three different two- 
body potential models. Figures IS.2I (a) and (b) show 
the results for the potential models v^(X = A|) and 
v a( A = A io) I Ec l- (E2J|], respectively, while Fig. E2 (c) 
shows the results obtained for v sc h{D = D$) [Eq. f|S . 1 [) ] . 
The most striking aspect of these figures is that U V (R) 
and W V (R) are substantially different, meaning that the 
nonadiabatic couplings P vv i (R) and Q vv < (R) [Eqs. (|S.13|) 
and (|S. 14[) ] are important near R = r v( jw- Therefore, 
it is clear that one needs to go beyond a single chan- 
nel approximation in order to describe the three-body 
observables. It is worth noting that the nonadiabatic 
couplings originate from the hyperradial kinetic energy. 
Their large values are thus consistent with our physical 
picture in which the local kinetic energy increases and 
generates the repulsive barrier in our effective potentials. 
It is also worth mentioning that since the three-body 
repulsive barrier prevents particles from approaching to 
small distances, the question of whether or not the short 
range physics actually changes as a function of the exter- 
nal magnetic field (as in experiments in ultracold quan- 
tum gases) |13j can not be directly answered by observing 
features related to Efimov physics. 

The strong multichannel nature of the problem can be 
illustrated by comparing the results obtained from a sin- 
gle channel approximation to Eq. (IS.10[) . i.e., W uv i{R) = 
[y 7^ y'), with our solutions of the fully coupled sys- 
tem of equations. Figure IS. 31 (a) shows the three-body 
parameter k* [related to the energy of the lowest Efimov 
state through the relation k* = (mE/h 2 ) 1 ^ 2 } obtained for 
the model obtained in the single channel approxima- 
tion (open triangles) as well as our full numerical results 
(open circles). The disagreement between these quanti- 
ties increases with the number of s-wave bound states n, 
meaning that the physics controlling the results becomes 
more and more multichannel in nature. Nevertheless, we 
find that by imposing a simple change in the adiabatic 
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(a) v a x (X = A*) 
4 




(b) »5(A = Aio) 
— 4 




(c) v sch (D = D* 6 ) 
„ 4 



ST 




FIG. S.2: Comparison of U V (R) (green dashed curves) and 
W V (R) (red solid curves) demonstrating the importance of 
nonadiabatic effects introduced by Q„ U (R) in Eq. (|S.lf [) . (a) 
and (b) show the results for the potential models v1(\ — A|) 
and v%(\ = Ai ) [Eq. (|S.2[) ]. respectively, and in (c) we show 
the results obtained for v sc h(D = Dg) [Eq. (|S.l[l ]. (a) and 
(b) also show the effect of the diabatization scheme used in 
order to prepare some of our figures in the main text (dash- 
dotted curves). The goal of the diabatization is to eliminate 
the sharp features resulting from Q„„(i?) in Eq. (|S . 1 f I) which 
should not contribute substantially to the three-body observ- 
ables. The case shown in (c), however, does not allow us to 
easily trace the diabatic version of the potentials W v . In this 
case, however, the "less" sharp features have a larger contribu- 
tion due to the crossing with a three-body channel describing 
a collision between a <?-wave molecular state and a free atom, 
giving rise to the anomalous n — 6 point for the v sc h model 
in Fig. 4 of the main text. Although such cases are relatively 
infrequent in our calculations (and occur mostly for the f sc h 
model), they do illustrate nommiversal effects that can affect 
the three-body parameter. Nevertheless, the three-body ob- 
servables obtained for cases like the one shown in (c) are still 
within the 15% range we claimed for the universality of the 
three-body parameter. 



potentials near the barrier — to make the barrier more 
repulsive [see Fig. IS. 31 (b)] — the single channel approxi- 
mation for [filled circles in Fig. IS. 31 (a)] reproduces the 
full numerical calculations much better. This agreement 
indicates that most of the nonadiabaticity of the problem 
is related to the exact shape of the barrier and that, to 
some extent, the effect of the nonadiabatic couplings is 
to make the effective potential W v more repulsive. For 
these reasons and, of course, the universality of our full 
calculations (see for instance Fig. 4 of our main text), 
we believe that the short-range barrier in the three-body 
effective potentials indeed offers a physically valid expla- 
nation of the universality of the three-body parameter. 
We emphasize, though, that Fig. 4 in the main text only 
includes the results from our essentially exact solutions 
of the full calculations. The single channel results dis- 
cussed here are intended only as support of our physical 
interpretation. 

It is within this context that we analyzed the model 
proposed in Ref. [l4|. In Ref. the three-body ef- 

fective potential important for Efimov physics was es- 
timated by considering the different aspects controlling 
the physics at small and large distances. At distances 
comparable to R = r v( jwj it was assumed that the ef- 
fective three-body potential is dominated by the contri- 
butions from equilateral triangle geometries and included 
only on two-body interactions. Under these assumptions, 
^12 = ?"23 = 7"3i = t giving R — 3 1 / 4 r (note that our 
definition for R differs from that used in Ref. [Hj]), the 
effective potential can be written as 



V m (R) 



-C 6 /rf 2 



Ce/4 3 - C e /r e 31 



-3C 6 /r 6 = -3 x 3 3/2 C 6 /R 6 - (S.15) 



This potential is expected to be valid for distances R < 
A = 47rr(l/4)- 2 3 1 / 4 3 3 / 8 r vdw « 1.9r vdw [2|. With our 
method, however, we have the means to check the valid- 
ity of Eq. (|S.15[) by comparing it with our numerically 
calculated potentials. Figure IS. 41 shows the three-body 
potentials obtained using the v^(X — A 10 ) model support- 
ing a total of 100 two-body bound states. The potential 
from Eq. (|S.15[) is the black solid line passing near the 
series of avoided crossings and might be loosely viewed 
as diabatically connecting the fully numerical potentials. 
This relation is reasonable given that this sequence of 
avoided crossings has been shown in Ref. [l5[ to be re- 
lated to the transition of the system from an equilateral 
triangle geometry to other geometries. This figure there- 
fore suggests that approximating the short range physics 
by Eq. (|S.15[) is not wholly unphysical, but its validity 
depends on a strong assumption of diabaticity through a 
large number of avoided crossings and is thus probably 
not an approximation satisfactory for a more quantitative 
analysis. 

The potential in Eq. (|S.15[) . however, was not actually 
used in the calculations in Ref. 14j . Instead, it was used 
to motivate a claim that three-body quantum reflection 
plays an important role, allowing the short range behav- 
ior to be replaced by a repulsive potential for distances 
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FIG. S.3: (a) This figure compares the energies (as char- 
acterized by the three-body parameter ««) obtained from a 
single channel approximation with our full calculations. The 
three-body parameter k, is shown for the 1>J model in the 
single channel approximation (open triangles) as well as for 
our full numerical results (open circles). The single channel 
approximation can be improved by imposing a simple change 
in the adiabatic potentials near the barrier, as is shown in (b). 
There we smoothly connect the potential for v% (red solid line) 
to the barrier obtained for «vdw (black solid line), resulting in 
the potential labeled by v% (green solid line). This new poten- 
tial is actually more repulsive and has energies [filled circles 
in (a)] that are much closer to our full numerical calculations. 



R < A rj 1.9r v( jw- It is interestin g to note that the value 
R s» 1.9r v dw obtained from Ref. jl4[ for the position of 
the hard wall is quite close to the hyperradius where our 
potentials exhibit the universal barrier (see Fig. IS. 21 for 
instance), indicating that A might have some physical 
meaning. It is worth mentioning, however, that the bar- 
rier we observed in our calculations is model independent, 
i.e., it doesn't rely on the particular model used for the 
two-body interaction. For distances R > A, the model in 
Ref. [lij assumed the three-body effective potentials to 
be given by the universal Efimov formula, 



V E (R) 



s 2 



1/4 



2[iR 2 



(S.16) 



It is well known [27|, however, that this potential is valid 




0.6 0.8 1 

Y?/r vdw 

FIG. S.4: This figure shows the three-body potentials 
obtained using the Vx(X = A*o) model supporting a total 
of 100 bound states. Roughly speaking, the potential of 
Eq. (pT5| pj (black solid line) can be seen as a diabatic 
potential since it passes near one of the series of avoided cross- 
ings. 



for r v( jw <SC \a\, and R = 1.9r v dw is certainly out of 
this range. In fact, one can see in Fig. 1 (b) of the main 
text, that the use of the Efimov potential [Eq. (|S.16[) ] 
for R < 10r v dw is a crude approximation. Nevertheless, 
using this model, Ref. [l4| obtains a j, » — 9.48r v dw> a 
value consistent with experiments (Iq - l26l | as well as with 
our present calculations for a^ h . However, extending the 
model of Ref. [l4| to the limit \a\ — oo, we obtained 
k* ~ 0.037/Vvdw- Our result, by way of contrast, using 
«vdw [ E q- dS3|] is = 0.226(2)/r vd w- We also have 
tested the effects of finite a corrections on the model of 
Ref. [lH by replacing Eq. f)S. 16[) with the three-body po- 
tential obtained with a zero-range model of the two-body 
interactions [28| . These corrections are particularly im- 
portant near R = \a\. With this modification, the model 
of Ref. " 



14] 



leads to a^ h w — 39.96r v dw- For these rea- 
sons, we believe that this model's agreement with our 
results and experimental data is fortuitous. 

While multichannel couplings are needed to quantita- 
tively describe this system, it is possible to construct 
an effective hyperradial potential curve that correctly 
describes the behavior of three atoms in the universal 
regime. Such a potential curve could be useful for sim- 
plified future studies. The approximate form obtained 
from the present study is: 



^f^W-(R) « - 

2.334 

(i?/r vd w) 3 
44.52 



4 



- V 4 



(-R/r v dw) 

1.348 

(i?/r vdW ) 4 
4.0 x 10 4 



44.52 4.0 x 10" 

~ (i?Avdw) 5 + (#Avdw) 16 ' 



(S.17) 
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A/r vd 




Here \i is the three-body reduced mass and r v dw is the 
two-body van der Waals length. For comparison the spec- 
ulative potential curve proposed by Chin [14] is shown, 
which does not resemble the present result even quali- 
tatively at small distances, as it exhibits far too little 
attraction in the region R = 2 — 5r v( jw- 



FIG. S.5: Comparison between the effective potential ob- 
tained from Ref. [Ij] and the potential from Eq. (|S.17[) . 
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